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ABSTRACT 

We consider photon emission from a supercritical accretion disk in which photons 
in the inner regions are trapped and advected in towards the center. Such disks are 
beheved to be present in failed supernovae and gamma-ray bursts, and possibly also 
in ultra-luminous X-ray sources. We show that the luminosity from a supercritical 
accretion disk is proportional to the logarithm of the mass accretion rate when the 
vertical profile of the matter density is exponential. We obtain analytical results also 
for other density profiles, but the dependence of the luminosity on the accretion rate 
is no longer logarithmic. 
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1 INTRODUCTION 

There is wide consensus that active galactic nuclei (AGNs), 
gamma-ray bursts (GRBs), X-ray binaries, etc., are powered 
by accretion flows on to compact relativistic objects, most 
often a black hole (BH). A well-known model of such accre- 
tion flows is the standard-disk model proposed by Shakura 
& Sunyaev (1973), which has been used to explain a variety 
of observations of AGNs and X-ray binaries. 

Recently, it has been recognized that observations of 
many objects cannot be explained by the standard-disk 
model. For instance, this model cannot produce the very 
high temperature {T > 10^" K) and broad-band spectrum 
(extending from 10^ to > 10^* Hz) observed in the Galac- 
tic Center source Sgr A* and in X-ray binaries in the hard 
state. This has led to the idea of an advection-dominated ac- 
cretion flow (ADAF, Narayan & Yi 1995, Abramowicz et al. 
19950) , which is also called a radiatively inefficient accretion 
flow (RIAF). The RIAF model explains the low luminosity, 
high temperature and optically-thin emission of these sys- 
tems. 

The standard disk cannot also be applied to very high- 
luminosity accretion disks in which the mass-accretion rate 
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^ We note that Ichimaru (1977) proposed similar ideas 20 years 
earlier. 

^ For recent developments on the RIAF model of Sgr. A*, see 
Manmoto, Kusunose & Mineshige. (1997), Yuan, Cui & Narayan 
(2005) and references therein. 



M exceeds the critical mass-accretion rate Merit (= L^/c? , 
where Le is the Eddington luminosity). In this regime, the 
disk becomes optically very thick and as a result photons 
are partially trapped inside the accreting gas. This is the 
regime of interest to us. 

Photon-trapping is important at any radius where the 
accretion time scale is longer than the photon diffusion time 
scale from the disk midplane to the surface (Begelman 1978). 
The photons produced in this region of the disk are advected 
into the central BH and are unable to escape from the flow. 
Begelman considered only spherically symmetric accretion, 
but it is recognized that photon-trapping is important even 
for systems with a disk geometry. 

Photon-trapping has been included approximately in 
the so-called "slim disk model" proposed by Abramowicz, 
Igumenshchev & Lasota (1998), and in recent calculations 
by Szuszkiewicz et al. (1996) and Watarai et al. (2000), as 
well as in numerical simulations by Honma et al. (1991), 
Szuszkiewicz & Miller (1997), and Watarai & Mineshige 
(2003). However, these studies do not model the effect fully 
since they consider vertically integrated quantities in the 
disk and reduce the problem to a one-dimensional model. 
Full two-dimensional radiation-hydrodynamical simulations 
(2D-RHD) were done for the first time by one of the current 
authors in Ohsuga et al. (2005) 0. 



Eggum et al. (1988) performed the first 2D numerical simu- 
lations of supercritical accretion disks, and their work was later 
improved by Okuda (2002). However, the computation time in 
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A realistic model of supercritical accretion must include 
the following two effects: i) Even at disk radii inside the trap- 
ping radius, photons that are emitted near the surface of the 
disk can escape, though most of the photons emitted deeper 
inside are advected into the BH. ii) The vertical profiles of 
various quantities in the disk such as density, optical depth, 
etc., play a critical role in determining what fraction of the 
photons escape, and thereby influence the total luminosity 
of the disk. 

The first effect was pointed out and studied by Ohsuga 
et al. (2002) using simple analytical models and numerical 
simulations. In this paper, we include the second effect and 
study the combined influence of both effects for a broad class 
of accretion-disk models. We develop an analytical model 
and show that it agrees well with the numerical results of 
Ohsuga et al. (2005). 

The results obtained here may change the current un- 
derstanding of models such as the slim-disk and convection- 
dominated accretion flow model (CDAF, Stone, Pringle 
& Begelman 1999; Narayan, Igumenshchev & Abramowicz 
2000; Quataert & Gruzinov 2000; Igumenshchev, Abramow- 
icz & Narayan 2000), as well as the neutrino-dominated 
accretion flow (NDAF) model (Popham, Woosley & Fryer. 
1999; Narayan, Piran & Kumar 2001; Kohri & Mineshige 
2002; Di Matteo, Perna & Narayan 2002; Kohri, Narayan & 
Piran 2005 Gu, Liu & Lu, 2006) which is believed to play 
an important role in the central engine of GRBs (Narayan, 
Paczynski & Piran 1992). Q 

In the following we use R to represent the cylindrical 
radius of a point with respect to the central mass and z 
to represent its vertical distance from the equatorial plane 
of the disk. We also assume that the temperature T in the 
disk is lower than the electron mass {kT < meC?), so that 
the electrons are nonrelativistic. 



2 STANDARD ONE-DIMENSIONAL 
APPROXIMATIONS 

In this section we introduce the standard one-dimensional 
(ID) treatment of accretion disks, following the approach 
described in Ohsuga et al (2002) . In this approach, we inte- 
grate along the 2-axis and omit the ^-dependences of various 
astrophysical quantities. For instance, we write the optical 
depth from the mid-plane of the disk [z — 0) to the surface 



(^Z — 2max) as 



(1) 

(2) 



where a^e is the scattering cross section of a photon off an 
electron. Although a^e depends in general on the photon en- 
ergy E~i, we limit ourselves to Thomson scattering, for which 
the cross-section is independent of energy. This is reasonable 



these studies was short and the numerical models did not reach 
steady state. Ohsuga et al. (2005) were the first to successfully 
compute models approaching steady state. 

^ The relativistic magnetohydrodynamic (MHD) jet model is an- 
other attractive model of energetic accretion sources (e.g., McK- 
inney 2005; 2006 and references therein). This model is less likely 
to be affected by our work. 



for the low energy photons <^ 0{MeV) of interest to us. 
H is the disk half-thickness and is the electron number 
density in the disk mid-plane. 

The diffusion time scale for a photon to escape from the 
disk mid-plane is given by 



(3) 



where A'^scatt is the number of scatterings, A — l/((TTWe) is 
the mean free path, and c is the speed of light. Since Afgcatt 
is given by 

A^scatt « 3 r^ (4) 

(assuming a random walk in three-dimensions, see Ap- 
pendix |B]|, we write 

tdiff « 3tH/c. (5) 
The accretion time scale is 



t --^ 

f-acc — 1 
Vr 



(6) 



where the radial velocity Vr is related to the mass accretion 
rate M and the density p by 



M 



AtvR pdz 

M 
4nRpH ■ 



Thus 



2r R^ 
rh c Rn 



(7) 
(8) 

(9) 



where Rg = 2GM/(? is the Schwarzschild radius, and we 
have set p ~ nipUe and introduced a dimensionless mass 
accretion rate rh, 

rh = M/ Merit- (10) 

The critical mass accretion rate Merit is defined in terms 
of the Eddington luminosity Le as follows. Merit = L-e/c? , 
where Le = '^-KcGrUpM / aT ^ 1.3 x 10^** erg sec~^ (M/Mq) 
for a black hole of mass M. 

Photons can escape from the disk only if the condition 
tditi < ta.cc is satisfied. Therefore, the disk radiates freely 
only from radii larger than a certain limit, 

R > -Rtrap, (11) 

where the trapping radius i?trap is given by 



J?trap — -^hinrhRg . 



and 

hin = 



H 

"r 



(12) 



(13) 



flS: -Rtrap 

assuming that, when R < -Rtrap and radiation is trapped, 
the disk is geometrically thick with kin ~ 1 and that h 
here is independent of R. As we show in Appendix [X] h is 
indeed constant (« 0.5) for R ^ -Rtrap and decreases oc 1/R 
for -Rtrap < -R. This vahdates Eq. (fT^ . Note that -Rtrap is 
proportional to rh. 

The energy released in the disk per unit time by viscous 
dissipation is given by 

/•flout 

L = 2 / 2nRQ^i4R)dR, (14) 
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where the viscous heating rate is 



^ , 3 GM ,v 



(15) 



Here _Rin and -Rout are the radii of the inner and outer edge 
of the disk. However, not all the released energy is radiated 
because of photon-trapping. We discuss two cases below, 
depending on whether J?trap is larger or smaller than _Rin. 
For readers' convenience we introduce a dimensionless inner 
radius normalized by the Schwarzschild radius. 



(16) 



LlD,Rt,„.p<R + Lin,R<Rtra 

1 



2/lin 

1 

2h~ 



Le 



l + 2hi 



l + 2..1n(^) 

3/lin 



Rt) 



Rin 
27/lin 



-1/2' 



-1/2 



(20) 



Note that the luminosity is reduced relative to the result in 
Eq. (|17|l for the untrapped case, and increases only logarith- 
mically with rh. This explains why disks with supercritical 
accretion (m 3> 1, "slim" disks) are radiatively inefhcient. 



2.1 -Rtrap < Rin 



3 TWO-DIMENSIONAL TREATMENT 



When -Rtrap < -Rin, there is no trapped region and we esti- 
mate the total luminosity to be 



L 



lD,tot,.Rtrap<-Ri„ 



2 / 2TvRQ^is{R)dR 

1 -Rtrap 
6/lin -Rin 

Le . 
m, 

4rin 



(17) 



which is linearly proportional to rh. This is a well known re- 
sult for the standard thin accretion disk model in Newtonian 
gravity. 

2.2 -Rin < -Rtrap 

When -Rin < -Rtrap, only the region of the disk outside -Rtrap 
radiates freely and the luminosity from this region of the 
disk is 

iiD,.Rt„p<fl = 2 / 27ri?Qvis(i?)d-R 
J Rti„p 



2/lin 



2/lin 



1 - 



1 - 



2 / -Rtrap \ 

3 V i?in / 



8rin 
27/ii„ 



-1/2- 



-1/2 



(18) 



Within the one-dimensional approximation being discussed 
in this section, it is unclear how much luminosity is emitted 
from the region of the disk -Rin < R < -Rtrap- One extreme 
assumption is that there is no luminosity at all from this 
region, as in Ohsuga et al. 2002. Alternatively, one might 
assume that the disk flux here is limited to the local Ed- 
dington flux F « -Fe(-R) = LE/{4:nR^) = GMnipc/iR^aT). 
The luminosity from -Rin < R < -Rtrap of the disk is then 
given by 

I'iD,.R<.Rt„p = 2 / 2nRFEM{R)dR 



/3^in . \ 

m 

V 2rin J 



Le In I 



(19) 



If we include the term given in Eq. (|19[) . keeping in mind 
that it is uncertain, then the "total luminosity" becomes 



So far we have discussed a ID approximation and assumed 
somewhat arbitrarily that the region -R < -Rtrap either does 
not radiate at all or radiates at the local Eddington rate. 
In this section we carry out a more careful two-dimensional 
analysis. We focus on the region -Rin ^ R < -Rtrap discussed 
in Subsection 12.21 but now we carefully consider the two- 
dimensional distribution of quantities. 

In the following, we consider most of the quantities 
in the disk to depend on both R and z (the vertical 
height). However, we assume that the thin-disk approxima- 
tion is approximately valid in the sense that we omit the 
z-dependence when considering hydrostatic balance in the 
^-direction (cs — QkH , with the sound speed eg and the 
Keplerian angular velocity fix — ^ GMjR?) and angular 
momentum conservation (y J pdz — Mil — -Rin/-R)/(37r), 
with the kinetic viscosity u). Then the z-dependent optical 
depth at a given -R is given by 



T[z) — I dz a^ene{z) 
ctt / dz ne{z), 



(21) 
(22) 



where ne{z) is the electron number density at height z, and 
Zmax corresponds to the upper surface of the disk. We see 
that the expression for r given in Eq. ([T| in the ID approx- 
imation corresponds to r(0) in Eq. (|2ip . 

We now introduce a detailed model for the vertical 
electron density profile. As seen in Fig. [1] we can flt the 
the density profile found in the numerical simulation in 
Ohsuga et al. (2005) quite well with the following power- 
law form, 



-5^ 



q-l 



(g > 1, and ^ z qH), (23) 



where ne(0) and the power-law index q may depend on the 
radius -R. In the limit as g ^ oo, this model takes the form 
of an exponential 



Tie (z) = ne (0) exp ( - — 



(-J) (OsS^sSoo). 



(24) 



JlD,tot,Ri„<H,t,-a 



^ If on the other hand we have a constant 2-independent density, 
as considered by Ohsuga et al. (2002), it would correspond to 
q = I. The corresponding results can be recovered by setting 
g = 1 in all the expressions in the current paper. 



4 K. Kohri et al. 



sion for the ^-dependent diffusion time scale, 
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Figure 1. Variation of the matter density p with height 2 at a 
radius _R = 20-Rg. The solid line represents the result from the 
numerical simulation of Ohsuga et al. (2005). The data are fit well 
with either an exponential (dotted line) or a power-law (dashed) 
profile, with H/R = 0.4 and (for the power-law profile) q = 10. 
The numerical results in Ohsuga et al. (2005) can be fit with the 
same value of H/R = 0.4 for all radii up to ~ SOiJg. 



The optical depth from z to the surface is 

z 



r{z) = r(0) 1 



(25) 



except for g oo when it takes an exponential form. 

Because we adopt the numerical results corresponding 
to m > 100 in Ohsuga et al. (2005) , the radii R of order tens 
of Rg that we consider are smaller than -Rtrap ~ 10^ Rg and 
so we are in the radiation-trapped regime. Note that Fig. [T] 
gives H/R « 0.4 for the numerical simulation of Ohsuga et 
al. (2005), which is approximately consistent with the ap- 
proximation hin = 0.5 that we made in the previous section 
and in Appendix El 

The photon diffusion time scale from height z to the 
surface of the disk at Zmax is given by 

tdia{z) ~ A''scatt(z:, Zmax) Hz) / C, (26) 

is the number of scatterings in the 3D 



where Nscs.uiz, 
random walk, 

A^scattl^, Zmax) = 3 [t{z)]^ 



(27) 



and \{z) is the mean free path which is estimated using the 
properties of the disk at the starting position z of the photon 

1 



A(z) = 



O-T ne{z) 



The detailed derivation of the above relation is given in Ap- 
pendix[B] By using Eqs. (|27|l and (I28|) . we obtain an expres- 



tdiffiz) = St{z) \{z)/c 



3r(0) 



H 



r(0) 



(9 + 1)/'? 



(29) 



By requiring the diffusion time scale fdiff(z) to be 
shorter than the accretion time scale tacc given in Eq. (|9]), 
we find that photons can escape only from the region where 
the following condition is satisfied, 

9/(9+1) 



r(2) 
r(0) 



R 



Rtra 



(30) 



The photon-trapping radius i?trap is determined by Eq. (|12|l 
and is proportional m. A schematic picture indicating the 
region of the disk from which radiation can escape is shown 
in Fig. [2] We see that trapping is not determined by simply 
a critical radius 7?trap, but is described in terms of a two- 
dimensional surface ztra,^{R), where 

9/(9+1) 



R 



T(2trap) _ 

r(0) " \^i?trap 

Note that 2:trap is a function of R in this case given by 
1/(9+1) 



Zt-r. 



1 



R 



Rti. 



(31) 



(32) 



and its second-order derivative by R is always negative 
(f ztis.p / dR^ < for i? < -Rtrap because 

1/(9+1)' 



dztr 



dR 
and 

d? ZtT3.l 

dR^ 



qhin 



q + 2 



R 



Rti. 



hin 



Rti 



R 



Rti 



-9/(9+1) 



(33) 



(34) 



^ ~r -L ttrap \ i-trap 

for R < -Rtrap. The luminosity from the region Rin ^ R ^ 
Rtrap is then given by 

i2D,fl<flt„p = 2 / 2n RQ.isiR f ^y^^^ ' dR. (35) 



r(0) 



Here we assume that all the energy released at z ^ Ztrap is 
radiated, while the energy released at 2; < Ztrap is completely 
trapped. Also, because Qvis and t{z) are both proportional 
to J pdz, we express the vertical distribution of viscous dis- 
sipation directly by means of the factor r(2:trap)/T(0). 

Evaluating the above integral, we obtain the luminosity 
of the disk from the region R < -Rtrap: 



L 



R 



2D,.R<.Rtrap 
/■-Rtrap 

2 / 2nRQ^i,{R) 

J Ri„ \Wtrap 
1/(9 + 1) 



9/(9+1) 



dR 



^ r f Rtiap \ 

2hrn\R~) 



-Rtrap / R\i 



dx 



x(9+2)/{9+ 



)/{9+i) V a;i/2 



2{q + l) fR 



' + 3 



/ Rtrap \ 

V~R~"J 



-1/2 



+ 



(g + 1)' fRt 
q + 3 



f Jltrap \ 

V~R~"J 



1/(9+1) 
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Figure 2. Schematic representation of a supercritical accretion 
disk around a BH. The shaded region shows the photon-trapping 
region from which no photons can escape. We see that, even for 
R < -Rtrap, there exist regions near the surface of the disk from 
which photons can diffuse out (white region). 



2(g + l) /3/ii„m 



(36) 



' + 3 



2nn 



'''' , (g + if 



2nn 



1/(9+1) 



The dominant contribution comes from the second term, 
which is proportional to m^^^'''^^'^ for m ^ 1. The luminosity 
from the region 7?trap < R is, of course, the same as in 
the ID case (Eq.O, i.e., I/2D,flt™p<-R = iiD,flt„p<fl, since 
there is no photon trapping. Therefore, we have the following 
expression for the total radiated luminosity from the entire 
disk, 

-£'2D,tot(g) 

1 

4q 



-LeX 



3(g + 3) 
1 

' 2hi 
4q 



/ Rtva.p \ 
\ Rin J 



-1/2 



(g+1)' fRt 
g + 3 



^in 



1/(9+1) 



(37) 



3(q + 3) V 2n„ 



+ 3) V 



'''' , + 



V 2ri„ ) 



1/(9+1) 



In a similar fashion, we obtain the following result for 
the case of an exponential density profile (q oo). 



,tot((7 = oo) 
1 , 



2/lin 
1 

2h~ 



Le 



V Rin J 

32rin . 



-1/2 



-1/2 



hi 



+ ln 



Rti 



Rin 
(Shin^ 

\ 2n„ ' 



(38) 



To verify that this result is consistent with the limit q oo 
of Eq. (I36|l . we make use of the following useful relation. 



hm - 



2e + l 



In a:: 



(39) 



This shows that the dominant terms in Eq. (|36p do not di- 
verge as q — > oo but tend to a finite value that is proportional 
to ln(m) for m S> 1. 



4 COMPARISON WITH NUMERICAL 
SIMULATION 

Here we compare the analytical formulae derived above 
for the luminosity L of the disk with the results 
of multi-dimensional numerical simulations reported by 
Ohsuga et al. (2005). 

In Fig. [3] we plot the disk luminosity in Eddington 
units as a function of the dimensionless mass-accretion rate 
m. Analytical formulae are denoted by solid lines and cor- 
respond, in descending order, to the following cases: no 
photon-trapping (Ea. ll7|) . power-law density profile with q = 
1 and q = 10 (Eq. I37|l , exponential density profile (Eq. 1381 
i.e., q — » oo), and the standard ID model which ignores 
photons emitted from Rin < R < i?trap (Eq. I19p . Triangles, 
circles, and squares connected by long-dashed lines indicate 
the numerical results from Fig. 7 of Ohsuga et al. (2005) 
corresponding to metallicity Z — 0, IZq, and 10 Zq, re- 
spectively. The dotted line represents the results of the 
hydrodynamical simulations of the slim-disk model carried 
out by Watarai et al. (2000), which did not explicitly include 
the effects discussed here. In these calculations, we adopted 
Rin = 3i?g and h = 0.4. 

From Fig. [3] we see that the luminosity is larger than 
with the standard analytical ID model if we ignore radiation 
from the photon-trapped region. The increase is the result 
of the multi-dimensional nature of photon diffusion, which 
allows radiation to partially escape even from radii R < 

Rtrap • 

For reference, we have also plotted by a short-dashed 
line the approximate formula Eq. (|20|l for the disk luminos- 
ity in the ID model, which includes the photon luminosity 
emitted from Rin < R < -Rtrap by assuming a simple local 
Eddington fiux. Because of the logarithmic functional form 
of this approximation, this model resembles the exponen- 
tial density model in the 2D treatment. Correspondingly, 
we conclude that this model provides a reasonable ID ap- 
proximation. 

Ohsuga et al. (2005) pointed out that photons tend to 
have momenta pointed inwards because of advection by the 
accretion flow, which causes the radiation to spend a longer 
time in the disk before escaping relative to the simple dif- 
fusion estimate we have used in our model. However, we 
believe this introduces only a minor correction to our re- 
sults because we are focusing mainly on photons for which 

idiff < tacc. 

In addition, Ohsuga et al. (2002) obtained outflows 
driven by the super Eddington luminosity of their models. 
Although our analytical model does not include outflows, we 



^ The existence of metaUicity has two conflicting effects which in- 
fluence the disk structure, i) larger radiation coefficient induced 
by the sizable amount of metals helps to cool the gas more, which 
means that the gas can be easily accreted, and ii) the larger ab- 
sorption coefficient by such a metallicity pushes the gas outside, 
which means that it is difficult to be accreted. However there is 
an uncertainty of metallicity in the system. Therefore, Ohsuga et 
al. (2005) studied three cases of the metallicity in a reasonable 
range, i.e., Z/Zq=0, 1, and 10. The variation of the luminosity 
which comes from this uncertainty should be treated as an error 
in the numerical simulation. For details see Ohsuga et al. (2005). 
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Figure 3. Variation of the disk luminosity as a function of the 
mass-accretion rate. Analytical formulae are denoted by solid 
lines. Prom above, the lines correspond to the case with photon- 
trapping ignored, power-law density profiles with 9 = 1 (constant 
density), q = 10, and g — > oo (exponential profile), and a standard 
ID model in which photons are assumed to be fully trapped for 
-Rin < R < Rtmp ■ The triangles, circles, and squares connected by 
long-dashed lines indicate the numerical results shown in Fig. 7 
in Ohsuga et al. (2005), for metallicity Z = 0, Zq, and 10 Zq, 
respectively. The dotted line represents the case of the hydrody- 
namical slim-disk simulations of Watarai et al. (2001), without 
explicitly including 2D effects. The short-dashed line shows the 
profile of luminosity in an approximate analytical ID model in 
which we assume that the disk emits at the local Eddington rate 
at radii < R < i?trap. AH results correspond to = 3Rg 
and h = 0.4. 



believe our discussion is still applicable in a time-averaged 
mean sense. 



5 DISCUSSION AND CONCLUSION 

In this paper, we have studied photon emission from a super- 
critical accretion disk. In the usual one-dimensional treat- 
ment of disks, one integrates all quantities along the verti- 
cal (z) direction. This is reasonable for a thin accretion disk 
when the accretion rate is well below the Eddington rate. 
However, we show that, when the accretion rate is large, 
it is quite important to consider the z-dependence of fluid 
quantities and to detreat the vertical diffusion of photons 
carefully. If we omit the z-dependence of the diffusion pro- 
cess in the simplest version of the one-dimensional model, 
we tend to significantly underestimate the luminosity. On 
the other hand, if we omit the ^-dependences of the matter 
density and assume uniform density, we overestimate the 
luminosity. Only when we consider the z-dependence of the 



matter density and adopt the appropriate density profile 
such as a power-law or exponential form as a function of 
z, as seen in 2D radiation-hydrodynamical simulations, are 
we able to reproduce the luminosity found in the numerical 
simulations. 

Using the analytical approach developed here, we 
should be able to model how the peak energy of the photon 
spectrum decreases when the photons produced near the BH 
are trapped, as pointed out by Ohsuga et al. (2003). This 
will be discussed in a forthcoming paper (Kohri & Ohsuga, 
2007). 

When we simultaneously include radiative, advective, 
convective, and neutrino cooling processes in an analyti- 
cal model of the accretion disk, the effects discussed here 
will play a crucial role. In neutrino dominated accretion 
flows (NDAFs), for instance, the neutrino luminosity and 
annihilation rate could be modified substantially relative to 
the predictions of one-dimensional models. The dominant 
cooling process could be changed, and the entire picture of 
the accretion might be dramatically modified, e.g., for the 
produced energy through the v v annihilation (Di Matteo, 
Perna & Narayan 2002, Gu, Liu & Lu 2006, Chen & Be- 
loborodov), r-process nucleosynthesis (Surman, McLaugh- 
lin, Hix 2006), and so on. We plan to discuss these effects in 
forthcoming papers. 

Multi-dimensional hydrodynamic simulations of accre- 
tion disks including all of these cooling processes have not 
been done so far (except for the 2D numerical simulations 
of neutrino cooled disks by Lee & Ramirez- Ruiz, 2006). To 
clarify the role of multi-dimensionality on CDAFs, NDAFs, 
etc., it would be useful to compare the analytical estimates 
obtained in this paper with full multi-dimensional numerical 
simulations. It is hoped that these simulations will be done 
in the near future. 
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APPENDIX A: ANALYTICAL ESTIMATES OF 
DISK-HALF THICKNESS H 

In this section, we simply try to analytically estimate the 
i?-dependence on disk-half thickness H . The force balance 
along z-axis is generally expressed by 



F 



GM H 



(Al) 



with a photon flux of F. 

If ^trap < R, then it would be reasonable to assume 
that the flux is equal to the viscous heating rate, F ~ Qvis. 
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Then, we obtain 



H_ 

R 



3R^ 
4 R 



(A2) 



Therefore H/R is approximately proportional to oc 1/R for 

Rtrap < R- 

On the other hand, if J? < -Rtrap, we may assume that 
the flux would be approximately the order of the Eddington 
flux F « i^E = iLE/(47ri?^). The factor 1/2 in front of Fk 
is attached as a matter of convenience for the consistency. 
That comes from the viewpoint of continuities of astrophys- 
ical quantities although physics does not change at all by 
this artificial factor. From Eq. ([Aip we find that H/R is 
constant for R < -Rtrap, 

H 

"r 



(A3) 



Then we see that h = H/R is a smooth function of R, 



( 1/2 {R ^ iJtrap), 

h=l 

I ■2Rtia.p/R {Rtia.p < R), 

if we take -Rtrap = 3/2hinrhRg with hin = 1/2. 



(A4) 



APPENDIX B: DERIVATION OF THE 
Z-DEPENDENT /INDEPENDENT DIFFUSION 
TIME SCALE 

Here we discuss the diffusion time scale along with z-axis 
in 3D random walk processes when the scattering length 
depends on the position, i.e., in the case that we consider the 
2-dependent number density of electron, the cross section, 
and so on in the electron-photon scattering processes. 

In general, we can write the summed distance measured 
along the 3D random walk process from the stating point 
2 = zo as 



d{zo, 



E 







dz. 



(Bl) 



where -/Vscatt (zo, Zmax) is the number of the scattering from 
z = 20 to a maximum of z (= Zmax). The integral interval 
[zi^-i, Zi^ \ is determined by solving the following differential 
equation for the number of the scattering, 

dN 
dz 

When we solve it, we may assume that the interval between 
the scatterings is determined by the approximate relation, 



ene{z). 



(B2) 



,ne{z)dz 



(B3) 



where the meaning of dividing by three comes from the con- 
tribution only to 2-ajcis in the the 3D random walk. Of course 
the right hand side would not have to be one third if only 
it were the order of 0(1). For simplicity, here we have just 
took it one third. In addition, generally Zi^ must not be 
larger than Zi^-i. 

As we will see later, the diffusion tends to be proceed 
outward with increasing diffusion length or mean free path. 
Then the time spent in inner regions tends to be much longer 



than that in outer ones. That means that phenomena of the 
diffusion which started from 2 = 20 are mainly determined 
by the local physics at around 2 = 20. Then Eq. (|B3P might 
be rewritten as a definition of the 2-dependent mean-free 
path \{z) through 



by 



\{zi 



a^ene{z)dz ^ (7Tne(2i^_l) 



dz 



dz, 



3 (TT ne{Zi^-l) 



(B4) 



(B5) 



where we assumed that a^e ~ ctt- Although there would be 
a lot of ways to define the 2-dependent mean-free path, the 
definition in Eq. HB5|I would be relatively natural in the cur- 
rent context of the accretion disks because the phenomena 
are mainly locally determined. 

Then d(2o,2niax) in Eq. (IBip is rewritten as 



d(2o,2niax)= ^ l[zi^-i), 



with 



A(2iJ (2i, + l > 2jJ, 
— \{zi^) (Zi^ + l < Zi^). 



(B6) 



(B7) 



Next let us consider the averaged value of d{z) after 
sufficiently a lot of tries. 

-fVacattCzO.Zmax) 

(d(20,2nrax)) = ^ {£{z^^-l)), (B8) 

where (X) means the average of X after such trials. Here we 
can assume {£{zi^-i)) « 0. That is because we have assumed 
that the local physics at around 2 — Zi^-i determines the 
phenomena, and surely then this approximation would not 
be so bad. Then we see that the averaged value of d{zo, 2max) 
vanishes, 

(d(20,2nrax)) «0. (B9) 

On the other hand, however, the averaged value of the 
square of d{z) must be finite. 

(^d{zo, 2max)^) 

-fVscatt (ZO . Zmax ) -iVgcatt ( ZQ , Zmax ) 

-WscattCzo.Zmax) 

^ {£{z,^-,f) + w^^.-inz,^-!)) 

-WscattCzo.Zmax) 

« E (^(^«.-)^), 

where we approximated 



(BIO) 



J2 (^(2,.-l)£(2,,_l)) 



izT^jz \ izT^jz / 
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0, 



(BU) 



McKinney, J. C. 2005, arXiv:astro-ph/0506367 



because the step iz does not correlate with that of jz among 
the independent trials for {£{zi^-i)} = 0. Now we can ap- 
proximate (^X{zi^-i)^^ « A(zo)^ because the scatterings at 
the outer regions do not frequently occur and hardly con- 
tribute to the summation at all. 

Then from ifBlO)) . we find that the number of the 
scattering of the photon diffused from z to ^max, i-e., 
(d^(z,2max)) — (/^^"""' dz)'^ , IS represented by[3| 



TVs, 



t (-2 , 2niax) 



3 [aTUeiz)] 



dz 



f7T 



dz ne{z) 



= ir[zf 



(B12) 



where we have used the same logic in Eq. (|B4|) to transform 
the first line to the second. 

Here we get the expression of the z-dependent diffusion 
time scale, 



tdis{z) = A''scatt(2:, Zmax)A(2:)/c 
« iT{zf\{z)/c. 



(B13) 



Of course, if we omit the z-dependence, we immedi- 
ately get the z-independent diffusion time scale shown in 
Eqs. @ and ©. 

When we consider the power-law density profiles cx 
{1 — z/qHY"^ (or the exponential one as their large-g limit), 
surely tdiff(z) decreases rapidly as oc (l — z/qH)'''^^ as a func- 
tion of z. This means that the time spent in the inner regions 
is much longer than that in the outer ones, that validates 
our assumption in the current context in the accretion disks. 
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